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Analysis  of  Turbulent  near  Hakes 


A.  Prabhu  and  V.C.  Patel 


Abstract 

This  report  Is  concerned  with  the  analytical  treatment  of  the  flow  in  the 

i 

near  wake  behind  sharp  trailing  edges  with  a  fully-developed  zero-pressure 

gradient  turbulent  boundary  layer  just  upstream  of  the  trailing  edge.  The 

three  cases  considered  are  the  flows  behind  two-dimensional ,  axisymmetrlc  and 

infinitely-yawed  sharp  trailing  edges.  The  principal  limits  that  govern  the 

1/2  1/2 

near  wake  seem  to  be  c^  -*•  0  and  c^  Rfi  «*,  where  cf  and  Rfi  are  the  skin 
friction  in  the  upstream  boundary  layer  and  Reynolds  number  based  on  shear 
layer  thickness,  respectively.  This  Is  the  same  limit  as  that  for  the  local 
asymptotic  analysis  of  a  fully-developed  turbulent  boundary  layer  under  zero 
pressure  gradient  at  large  Reynolds  numbers.  These  limiting  conditions  lead 
to  the  governing  equations  for  the  inner  and  outer  near-wake  and  a  logarithmic 
variation  of  the  velocity  In  the  overlapping  region.  It  Is  argued  that  the 
same  limiting  conditions  also  lead  to  a  logarithmic  variation  of  the 
centerline  velocity  In  the  near  wake.  Similarity  solutions  have  been  obtained 
for  the  governing  equations  of  the  Inner  and  outer  near-wakes.  Limited 
comparisons  with  experiments  for  a  two-dimensional  flow  show  that  the 
variation  of  the  scales  and  profiles  predicted  agree  reasonably  well  with 
experiments. 


1.  Introduction 

Close  to  the  trailing  edge  of  an  arbitrary  body  the  flow  Is  usually 
dominated  by  separation  and  reattachment  or  trailing  edge  singularities  and 
this  complex  flow  Is  not  clearly  understood.  In  the  case  of  laminar  flow  on  a 
flat  plate  the  boundary  layer  on  the  plate  develops  Into  Goldstein's  near  wake 
beyond  the  trailing  edge  through  a  transition  In  a  small  region,  known  as  the 
triple-deck  region  (Stewartson,  1965)  around  the  trailing  edge.  This  triple¬ 
deck  region,  whose  order  of  magnitude  In  Reynolds  number  Is  higher  than  that 
of  the  boundary  layer  (Imbedded  Inside  the  boundary  layer)  arises  In  order  to 
avoid  the  singularity  of  the  boundary- layer  equations  at  the  trailing  edge. 
In  the  case  where  the  boundary  layer  Is  turbulent  the  extent  of  the  triple¬ 
deck  region  becomes  very  small  as  the  Reynolds  number  Is  necessarily  large  and 
therefore,  probably.  Is  of  little  consequence  for  the  development  of  the  near 
wake  beyond  It.  Evaluation  of  experimental  data  In  the  near  wakes  behind 
round  and  two-dimensional  bodies  have  been  recently  made  by  Patel  (1981a, 
1981b).  The  present  report  Is  concerned  with  the  analytical  treatment  of  near 
wakes  beyond  the  triple-deck  region  of  sharp-edged  bodies. 

If  the  oncoming  turbulent  boundary  layer  Is  fully  developed.  It  has  two 
distinct  layers  or  regions  governed  by  two  different  length  scales.  Since  the 
turbulence  structure  Is  known  to  retain  memory  It  Is  reasonable  to  expect  the 
flow  In  the  near  wake  also  to  retain  the  two-layered  structure.  Once  this 
two-layer  structure  (an  Inner  and  an  outer  layer)  similar  to  that  of  a 
turbulent  boundary  layer  Is  assumed  to  exist,  the  velocity  profile  In  the 
overlapping  or  matching  region  must  vary  logarithmically  as  In  a  turbulent 
boundary  layer.  The  logarithmic  velocity  variation  will  then  be  the  boundary 
condition  for  describing  the  Inner  and  outer  wake  development.  Further,  the 
changes  along  the  streamwlse  direction  x  near  the  trailing  edge  In  the  Inner 
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wake  are  largely  due  to  the  mixing  produced  by  eddies  whose  characteristics 
are  determined  by  the  viscous  scale  of  the  boundary  layer  near  the  wall. 
Therefore,  the  relevant  length  scale  for  development  In  the  x  direction  for 
the  Inner  wake  will  also  be  the  viscous  scale  of  the  upstream  boundary 
layer.  These  arguments  have  been  Incorporated  Into  the  study  of  two- 
dimensional  near  wakes  by  Alber  (1980),  Andreopoulos  and  Bradshaw  (1980), 
Ramaprlan  et  al.  (1981),  and  Robinson  (1967).  They  have  all  shown  that  the 
centerline  velocity  varies  logarithmically  for  large  x  In  the  near  wake. 
Alber  made  further  similarity  analysis  and  obtained  detailed  velocity  profiles 
In  the  Inner  wake  which  he  compared  with  the  velocity  profiles  measured  by 
Chevray  and  Kovasznay  (1969).  He  also  considered  a  region  of  development  very 
close  to  the  trailing  edge  similar  to  the  Goldstein  wake,  where  the  laminar 
sublayer  of  the  upstream  boundary  layer  Is  consumed.  There  exists  no  analysis 
for  the  outer  wake.  Also,  It  Is  not  clear  whether  similar  analyses  can  be 
made  for  near  wakes  behind  Infinite-swept  sharp  trailing  edges  and  pointed 
axl symmetric  bodies. 

In  this  report,  Goldstein's  near-wake  solution,  and  the  Inner  and  outer 
wake  analysis,  are  considered  for  the  three  cases:  two-dimensional.  Infinite- 
swept  and  axlsymmetrlc  wakes.  Goldstein's  1/3  power-law  for  the  variation  of 
the  centerline  velocity  seems  to  be  valid  for  all  the  cases  In  the  laminar 
region  very  near  the  trailing  edge.  Similarity  solutions  for  the  Inner  and 
outer  wake  which  have  logarithmic  velocity  variation  at  the  match  boundaries 
have  also  been  obtained.  Existence  of  similarity  solutions  In  the  Inner  layer 
requires  that  the  centerline  velocity  should  vary  logarithmically  for  large  x 
In  all  three  cases.  This  variation  Is  In  support  of  a  more  general 
dimensional  argument,  similar  to  that  of  Millikan  (1938),  which  leads  to  a 
logarithmic  velocity  variation  In  the  overlap  region  between  the  near  wake  and 


4 


the  far  wake.  Outer  layer  similarity  requires  that  the  square  of  the  total 
wake  width  vary  logarithmically  with  x,  and  this  Is  supported  by  the  recent 
measurements  of  Sastry  (1981). 


2.  The  Model  and  the  Equations 

The  various  flow  regions  In  the  near  wake  are  shown  In  figure  1.  The 
nomenclature  "near  wake"  here  signifies  the  region  beyond  the  trailing  edge 
over  which  the  memory  of  the  viscous  length  scale  of  the  boundary  layer  Is 
lost  or  where  the  log  region  of  the  boundary  layer  Is  consumed.  It  Is  not 
clear  at  this  stage  whether  the  asymptotic  or  the  far  wake  Is  far  beyond  the 
end  of  the  near  wake  or  matches  with  It  In  an  overlapping  region.  The  Inner 
layer  Itself  has  a  laminar  region,  very  near  the  trailing  edge,  where  the 
laminar  sublayer  of  the  upstream  boundary  layer  is  consumed  due  to  the  mixing 
produced  by  the  eddies  of  the  viscous  scale  and  a  region  beyond  it  where  the 
flow  Is  turbulent  but  characterized  still  by  the  dynamics  of  the  eddies  which 
retain  the  memory  of  the  viscous  scale  of  the  upstream  boundary  layer.  The 
relevant  length  scale  for  both  streamwlse  and  normal  direction 
Is  v/ut,  where  v  Is  the  kinematic  viscosity  and  u^  the  friction  velocity  of 
the  boundary  layer  just  upstream  of  the  trailing  edge. 

The  coordinates  chosen  and  the  corresponding  velocity  components  are 

shown  In  figure  2  for  all  three  cases.  The  governing  equations  for  the  outer 

and  Inner  flows  can  be  obtained  by  making  a  local  asymptotic  analysis,  similar 

to  that  made  for  a  fully-developed  turbulent  boundary  layer  (see,  for  example, 

Yajnlk  (1970),  and  Mellor  (1972)).  The  limiting  conditions  that  lead  to  the 

equations  given  below  are  the  same  as  that  for  the  turbulent  boundary  layer. 

They  are  u  /II  ♦  0  and  6u  /v  ♦  «.  where  6  Is  the  total  wake  thickness.  In 

T  u."u  2 

the  outer  layer  the  variables-^ — ,  t/u*,  y/5  and  x/L  are  held  fixed,  and  In 
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the  inner  layer  the  variables  u/ut,  t/ut»xut/v,  and  yuT/v  are  fixed  while  the 
above  limits  are  taken.  Note  that  for  the  boundary  layer  asymptotic  analysis 
(made  by  Yajnik,  1970,  Mel  lor,  1972),  the  stretched  variable  for  the  x 
coordinate  in  the  inner  layer  is  x/L  which  is  different  from  the  xut/v  used 
for  the  wake  here.  With  these  inner  and  outer  variables  fixed  under  the 
limiting  condition  u  /U  ♦  0  the  following  equations  which  govern  the  Inner 

T  ® 

and  outer  wake  are  obtained. 


Inner  wake 


(la) 

(lb) 

(lc) 


J  =  1  only  for  axl symmetric  and  k  s  1  only  for  yawed  wake,  otherwise  J  and  k 
are  zero. 


(2a) 

(2b) 


Here  u<j  and  w<j  are  the  defect  velocities  defined  by  u<j  *  Uw-u  and  w<j 
*  (W^-w),  and  t1  and  t2  are  the  shear  stresses  in  the  xy  and  yz  planes. 
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respectively.  Matching  between  these  two  layer  leads  to  the  boundary 
condition  In  the  overlapping  region  as 


where  u 

T 


u/uT  *  Ajln (-“*■)  +  Bj,  w/wT  =  Agin  (^JL)  +  Bg.  (2c) 

and  w  are  the  friction  velocities  In  the  x  and  z  direction. 

T 


3.  Laminar  Hear  Make  Region 

In  this  region  the  laminar  sublayer  of  the  upstream  boundary  layer 
develops  like  a  laminar  wake  before  It  is  consumed  by  the  mixing  due  to  the 
convected  small  scale  eddies  of  the  turbulent  boundary  layer  on  either  side  of 
the  plate.  The  variable  In  equations  (1)  are  normalized  appropriately  using 
the  velocity  scale  ut  and  the  length  scale  v/ut.  Superscript  +  denotes  these 
normalized  variables.  The  shear  stresses  Tj  and  Tg  are  given.  In  this  region, 
by  the  viscous  stresses  p3u/3y  and  u3w/3y,  respectively.  The  equations  which 
govern  this  region  can  then  be  written  as 


. »  +  +J »  , ,  +  +J» 

3  (u  y  )  3  (v  y  ) 


3x 


3y 


o  , 


and 


3x  3y  y  3y  3y 

k(u+*V 

3x  3y  3y  £ 


0  , 


(3a) 

(3b) 

(3c) 


At  the  edge  of  the  laminar  sublayer  u+  and  w+  vary  linearly  with  y+.  This  Is 
the  boundary  condition  for  large  y+  used  for  obtaining  solutions  of  equations 
(3)  In  the  laminar  region  of  the  near  wake.  At  the  center-line 
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~r  =  0.  In  view  of  the  linearly  Increasing  velocity  far  away  from  the 

3y 

centerline  Goldstein's  (1930)  asymptotic  expansion  for  the  laminar  wake  for 
small  distances  from  the  trailing  edge  seem  to  be  quite  appropriate  for  the 
solution  of  these  equations.  The  first  term  in  the  Goldstein  power  series 
expansion  for  the  wake  flow  close  to  the  trailing  edge  constitutes  the 
similarity  solutions  for  these  equations.  The  variables  for  the  similarity 
solution  are  given  by 

u+  *}x*  1/3  f'(n), 

«*  •  »  j  g(n).  (4) 

A  1  +  //  +\l/3 

and  n  ■  f  y  /(x  )  , 

where  A  is  given  by  the  relation  ~r~  =  which  is  satisfied  in  the  laminar 

T  + 

sublayer  of  the  boundary  layer.  If  v  is  obtained  from  the  continuity 
equation  (2a)  and  the  above  variables  from  (4)  are  substituted  into  equation 
(3),  the  following  equations  result 

nJf'2-  (2+J)  f"  /  nJf'  dn  -  (nJfM)'  =  0  ,  (5a) 

o 

k[(f'g  -  2f g * )  -  g"]  =  0  .  (5b) 

Here,  and  later,  the  primes  represent  differentiation  with  respect  to  the 
appropriate  variables.  For  a  two-dimensional  flow  (5a)  reduces  to  the 
equation 
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f'2  -  2ff "  -  f"  =0.  (6) 

If  f*  =  g  then  equation  (5b)  can  also  be  reduced  to  (6).  Equation  (6)  is  the 
same  as  that  obtained  by  Goldstein  (1930)  for  a  two-dimensional  flow.  For 
axisymmetric  flow  the  similarity  equation  is  slightly  different  and  is,  from 
(5a), 

nf'2  -  3f "  /  nf'dn  -  (nf")'  =  0  .  (7) 

o 

Equations  (5a),  (6)  and  (7)  have  linearly  increasing  solutions  for  large  n  as 
the  most  dominating  term  in  an  asymptotic  expansion  for  large  n.  The  boundary 
condition  for  large  n.  In  view  of  the  definition  (4),  Is  given  by 

f  -  9n  -  g  (8) 

and  the  inner  boundary  condition  is 

f(0)  =  fH(0)  =  g'(Q)  =  0.  (9) 

For  small  n  it  is  possible  to  develop  a  series  solution  for  these  equations. 
In  the  two-dimensional  and  yawed  wakes  it  is  the  same  as  that  given  by 
Goldstein,  namely 

f'(n)  -  g  =  aQ  +  a2n2/2  -  2a*  n4/4  +  ...  (10a) 


In  the  axisymmetric  case  the  series  solution  for  small  n  assumes  the  form 
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f  '(n)  *  aQ  +  -5-  n2  -  35- n3  +  •••  (!Ob) 

where  a0  is  a  free  parameter  to  be  determined  by  the  requirement 
that  f'(n)  and  g(n)  increase  like  9n  at  the  edge  of  this  layer.  Equations 
(5a)  and  (5b)  have  a  certain  invariant  characteristic  which  leads  to  the 
following  property.  If  F^(n)  and  F2(n)  are  solutions  of  (5a)  and  (5b)  then 
f  =  aFj(e)  and  g  =  aF2(c)  with  5  =  an  are  also  solutions  of  equation  (5a)  and 
(5b),  respectively.  In  view  of  this,  a  *  f'  (n)/F{(c)  •  Since,  for 
large  n,  f'(n)  varies  as  9n  and  Fj(c)  as  Fj'(»)c  it  is  easily  seen  that 

a3  =  9/Fj («).  (11) 

F£(«)  can  be  obtained  by  constructing  a  solution  starting 
with  Fx (0)  *  Fj (0)  =  0  and  FJ(0)  =  1.  Since  aQ  *  f'(0)  =  a2Fj(0),  equation 
(11)  leads  to  a  relation  between  a0  and  F£(»)  as 

a0  »  (9/(Fi'(»)}2/3.  (12) 

Once  aQ  is  found  from  (12),  the  centerline  velocity  is  obtained  from 

u+(0)  =  w+(0)/A  =  aQ  J  x+  ^ 3 .  (13) 

Equation  (13)  leads  to  the  centerline  velocity  variation  very  close  to  the 
trailing  edge,  and  equation  (10a)  and  (10b)  lead  to  the  velocity  profiles  near 
the  centerline. 
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4.  Turbulent  Inner  Near  Wake 

In  this  region  the  flow  is  completely  turbulent.  It  is  therefore  assumed 
that  the  viscous  stress  is  negligible  compared  to  the  turbulent  stress.  The 
latter  is  obtained  by  an  eddy-viscosity  assumption.  Then  if  equations  (1)  are 
normalized  as  In  section  3  the  following  governing  equations  are  obtained 


_/  +  +J\  .  .  +  +J » 

KuJL  )  +  3  (v  v  )  =  0 
+  +  * 
ax  ay 


+  iu_  +  +  3u_  _  JL  _9_  ,  +0  aux 
ax  ay  J  ay  ay 


-  o. 


(14a) 

(14b) 


where  vT  is  the  ratio  of  eddy  viscosity  to  molecular  viscosity.  The 
assumptions  on  the  form  of  the  eddy  viscosity  will  be  made  later. 
Since  6ut/v  tends  to  infinity  in  the  limit  for  which  equations  (1)  and  (2)  are 
derived,  it  is  reasonable  to  assume  that,  even  for  the  wake,  a  Reynolds  number 
based  on  ut  and  a  length  scale  h  describing  the  extent  of  the  inner  wake,  also 
tends  to  infinity.  The  limit  u/ut  ♦  0  corresponds  to  lnffiu^/v)  •  and 
hence  [In  h  uT/v]-1  ♦  0.  Here  u+  is  expanded  in  an  asymptotic  expansion  in 
terms  of  the  small  parameter  E  =  (In  h  uT/v)_1  in  the  form 


+ 

u 


S(x  )  =  f  (♦q  +  e  +  e  ^2  +  •••)» 


(15) 


where  k  is  the  Karman  constant,  4^,  etc.  are  functions 
of  x+  and  y+,  and  S(x+)  is  related  to  the  centerline  velocity.  S(x+)  tends  to 
become  large,  as  will  be  seen  later  because  eventually  at  the  end  of  the  near 
wake  It  must  be  of  the  same  order  as  the  velocity  ratio  u+  at  the  edge  of  the 
log  layer  of  the  boundary  layer.  It  is  also  argued  later  on  that  if  there 
exist  two  length  scales  for  describing  the  wake  development  in  the  x 


11 


f 

t 


direction,  the  viscous  length  in  the  near  wake  and  the  distance  from  the 
trailing  edge  in  the  far  wake,  then  matching  between  the  two  regions  requires 
that  the  centerline  velocity  should  become  large  compared  to  the  friction 
velocity  in  the  overlapping  region.  The  right-hand  side  of  (15)  then 
represents,  in  some  sense,  small  departures  from  the  centerline  velocity.  It 
is  further  assumed  here  that  and  ^  are  functions  of  the  similarity 
variable  $  =  y+/h+  only,  where  h+  is  huT/v.  Then  substituting  (15)  into  the 
continuity  equation  (14a)  gives 


+ 

V  = 


+ 


+  *  £ 

^  +  e*“  +..)*  + 

kCJ  0  0  1 

+  '  £ 

e2  /  2e$2  +  ..) 


(16) 


The  momentum  equation  (14b)  then  reduces  to 


The  above  equation  contains  terms  of  0(1)  and  0(e).  The  terms  that  are 
neglected  are  0f  0(e2)  or  smaller  (see  also  discussion  below).  In  order  to 
separate  the  terms  in  the  above  equation  the  order  of  magnitude  of  the 
coefficient  of  each  term  has  to  be  first  assessed.  The  order  of  magnitude  of 


X 


each  term  will  be  clear  after  certain  conditions  for  similarity  solutions  are 
evaluated.  If  and  ^  in  equation  (17)  are  to  be  functions  of  c  only,  all 
the  functions  of  x+  appearing  as  coefficients  in  equation  (17)  must  be 
constants.  In  articular 


h+SS'  *  c1# 


(18) 


ci  and  C2  can  be  set  equal  to  unity  without  loss  of  generality.  Then  the 
solution  of  (18)  leads  to 


_  1  ,  .  + 

S  =  y  1 n  h  +c3 . 


h+[ln  h++c3-l]  =  k2x  ,  (19) 


where  C3  is  a  constant  which  could  be  used  to  define  a  suitable  Initial  h+. 

These  solutions  reveal  the  character  of  the  centerline  velocity  variations  and 

the  limits  In  which  these  approximations  are  valid.  S  varies  like  e"1  and  h+' 

+  +'2 

like  e.  Under  these  conditions  h  S'  *  ke  and  h  =  k  e.  With  these  orders  of 

4.  +  • 

magnitude  for  h  S'  and  h  ,  It  is  clear  that  the  terms  that  are  neglected  in 

2 

equation  (17)  are  of  0(e  )  or  smaller.  An  equation  for  can  be  obtained  by 
collecting  terms  of  0(1),  as 


rr7*o>‘  *«;■ 1  * 0 

c  kh 


(20a) 


Collecting  coefficients  of  order  0(e)  leads  to  the  equation  for  4^,  as 


1  J  VT 


—  uv  -v  4;)'  ♦  -? — ?-  r 


kh 


2°  cJ  0 


(20b) 
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Equations  (20a)  and  (20b)  reduce  to  those  obtained  by  Alber  (1980)  for  two- 
dimensional  flow.  However,  the  present  derivation  Is  more  consistent  than 
that  made  by  Alber  (1980)  since  he  assumed  an  asymptotic  expansion  of  the 
form  <t>  «  $0  +  e(x)$^  after  obtaining  a  governing  equation  for  4  assuming 
that  t  Is  a  function  of  c  only  (which  Is  Inconsistent).  In  view  of  the 

particular  nature  of  e  In  this  problem  the  derivative  of  e  contributes  to 

2 

terms  of  the  order  0(e).  Hence  the  equation  for  remains  the  same  as  that 
obtained  by  Alber.  Equation  (20a)  can  be  further  simplified  to 


♦ 


I  I  I 
0 


(21) 


A  solution  for  the  above  equation  can  be  written  down  as 

1  5  kh+r  C  +  j  Ci  khu 

♦  i'(0  =  l-J-  exp  (-/  dc)3CA3+  /  kh  exp  f1  - - ^  dc2  <*,].  (22) 

C  Vj  T  0  0  T 

where  A3  Is  an  arbitrary  constant.  Since  *  0  at  c*  0,  equation  (22)  can  be 
further  simplified  to 

kh+  C  j  «  kh+^3 

♦  o'k)  “  ?2<exP  "  ;  *2  *  (23) 

V|C  O  C2  T 

Integration  of  (23)  once  gives 

C  kh+  i  Ci  kh+c, 

♦  i(e)  -  ^(0)  ♦  /  -“j  (/  (exp  -  f1  *3)  dcj  d^.  (24) 

o  J  o  c2  T 

For  a  yawed  wake,  If  an  expansion  for  w+  given  by 


A^sl+  r  (+o(b)  +  6  ’M5)  + 


+ 

w  = 


(25) 


Is  used.  It  can  be  easily  shown  that  ^(c)  Is  equal  to  ♦q(c)  given  by  equation 
(24)  with  J  =  0. 

No  assumption  has  been  made  so  far  concerning  the  form  of  the  eddy 
viscosity.  This  should  be  such  that  It  reflects  the  turbulent  structure  of 
the  Inner  wake  and  the  resulting  velocity  profiles  obtained  from  (24)  and  (15) 
should  be  such  that  the  matching  condition  (boundary  condition)  given  by  (26) 
Is  satisfied.  In  the  Inner  boundary  layer  Vy  varies  like  ky+.  Therefore  It 
is  reasonable  to  assume  In  the  wake,  at  least  In  the  logarithmic  region, 
that  vT/kh+  varies  like  t.  If  this  form  Is  assumed  throughout  the  Inner  region 
of  the  wake  It  leads  to  the  condition  that  the  eddy  viscosity  Is  zero  at  the 
centerline.  In  fact,  as  the  wake  develops  the  eddy  viscosity  at  the 
centerline  should  Increase  gradually  so  that  at  the  end  of  the  near  wake 
development  it  is  more  or  less  constant  across  the  entire  wake  width. 
Assumption  of  a  linearly  varying  eddy  viscosity  all  across  was  made  by 
Alber.  In  such  a  case,  even  though  the  velocity  profile  can  be  shown  to  vary 
logarithmically  for  large  c,  becomes  singular  as  ?  ♦  0.  In  order  to  avoid 
the  singular  nature  of  for  small  c  a  non-zero  eddy  viscosity  Is  required 
at  c  ■  0.  One  alternative  Is  to  assume  vy/kh+  ■  (c  +  a)  where  a  Is  of  the 
order  unity.  A  slightly  general  form  which  gives  the  required  behavior 
for  Vy  Is  given  by 

— 5 — j  *  1  -  p'/p  »  (26) 

vT/kh+ 

where  p  Is  chosen  such  that  p'/P  ♦  1  as  {  ♦  0  and  p'/P  ♦  0  rapidly  as  c  ♦  *. 
The  detailed  form  of  p  Is  still  unspecified.  This  form  Is  assumed  since  It 
facilitates  the  asymptotic  expansion  of  the  Integrals  In  equation  (24).  Under 
this  assumption  an  asymptotic  expansion  for  large  c  can  be  obtained  for  ^  as 


♦  '(C)  ■  ♦•(0)  *  ,  IkLTElP  Ifc  )dt  *  J  /  (/  !  pfT'.p  y  dF^FldCt 

o  C1  1  1  0  ?lCl°  o-  Plcl  *V  1  1 


‘  u-Pj/pip  (;  ‘  in  dtj, 

0  ?lcl  C1  1 


(27a) 


where 


■;  pdpTdF  • 


(27b) 


The  only  term  which  gives  a  logarithmic  behavior  for+'fc)  as  ?  ♦  »  is  the 
first  integral  on  the  right  hand  side  of  (27).  The  expansion  of  the  first 
integral  for  large  c  is  given  by 


/  (1;p,/p)pIdi:1  *  [(I-pVp)pI]  In  c  -  /  -Jr  (l-p'/pJpDln  c  *  (28) 

o  ?1  1  o  *1 


Since  p'  tends  to  zero  for  large  c  ,  p  tends  to  constant  for  large  c  ,  It  can 
easily  be  seen  from  (27b)  that  pi  tends  to  unity  for  large  c.  Therefore,  the 
first  term  on  the  right  side  of  (28)  is  log  c»  The  other  two  Integrals  in 
(27)  contribute  to  the  constant  term  In  the  expansion  of  ♦ '(c)  for  large  c. 
Thus,  the  expansion  for  ♦  '($)  for  large  c  Is  of  the  form 


♦  '(?)  *  j  In  t  +  constant  +  0(c_1) 


(29) 


It  Is  shown  here  that  the  solution  of  equation  (14b)  has  the  required 
logarithmic  behavior  which  matches  with  the  boundary  condition  In  the 
overlapping  region.  Irrespective  of  the  detailed  behavior  of  p.  The  velocity 
profile  thus  satisfies  both  the  conditions,  namely  zero  gradient  at  c  ■  0  and 
the  logarithmic  variation  for  large  c. 


Equations  (24)  and  (19)  and  (15)  give  the  complete  development  of  the 
velocity  distribution.  The  centerline  velocity  variation,  ucL  as  given  by 
(15)  and  (19)  show  that  for  large  x+,  u*L  can  be  expanded  In  an  asymptotic 
form  (only  the  dominant  terms  are  given  here)  as 

v  1  r  tn  x+  +  c4  (30) 

where  C4  is  a  constant  (which  depends  on  C3  and  4^(0)). 

This  logarithmic  behavior  for  large  x+  could  be  obtained  purely  from 
dimensional  arguments  In  exactly  the  same  way  as  the  logarithmic  behaviour  In 
y+  for  a  turbulent  boundary  layer  Is  obtained  (see  Millikan,  1938).  It  Is  a 
consequence  of  the  fact  that  In  the  near  wake  ucL  1s  a  function  of  x+  and  away 
from  the  near  wake  ucL/U<#  Is  a  function  of  x/L.  If  there  exists  an 
overlapping  region  where  these  two  functional  dependencies  match,  then  by  the 
same  arguments  used  by  Millikan,  It  can  be  shown  that  uc^  varies  as  given  In 
equation  (30)  with  the  constants  of  the  log  law  undetermined.  The  logarithmic 
variation  of  the  centerline  velocity  Is  therefore  a  requirement  of  the 
constraints  Imposed  by  the  two  different  similarity  coordinates  which 
Individually  govern  the  near  wake  and  the  wake  beyond  It.  The  fact  that  the 
dynamics  reflect  this  constraint  requirement  suggests  that  the  wake 
development  In  the  streamwlse  direction  also  can  possibly  be  described  by  * 
two-layer  hypothesis. 


5,  Turbulent  Outer  Wear  Mike 

Equations  (2)  that  govern  the  flow  In  this  region  are  In  fact  the  small - 
defect  wake  equations  (Townsend,  1956).  Denll  and  Landweber  (1979)  studied 


these  equations  in  connection  with  the  development  of  the  defect  layer  in  a 
thick  axisymmetric  turbulent  boundary  layer.  The  particular  cases  which  they 
disregarded  in  their  analysis  as  not  relevant  to  thick  axisymmetric  boundary 
layers  seem  to  be  the  most  relevant  for  the  present  case.  The  analysis 

carried  out  here  follows  the  lines  of  Denli  and  Landweber  (1979).  Because  of 
the  similarities  in  equations  (2a)  and  (2b)  the  general  features  of  the 
analysis  and  conclusions  on  the  development  of  u^  and  wj  are  similar.  Hence, 
only  equation  (2a)  will  be  treated  in  the  subsequent  analysis.  In  terms  of 
the  nondimensional  quantities  given  by 

_  _  _  _  2 

ud  =  Ud/U«*  x  =  X/L»  y  =  y/5i»  T  *  t/put,  (31a) 


(where  L  and  6^  are  the  length  scales  for  the  initial  wake  in  the  x  and  y 
directions)  the  mean  flow  equation  (2a)  reduces  to 


U 


«i  3ud 


r~  ~ 

U*  L  9  X 
T 


l  a 


(yJT ) 


-J  _ 

y  ay 


(31b) 


p  — 

Using  the  transformation  £  *  x,  equation  (31b)  can  be  written  as 


3ud  13  /-d— v 


u; 


where  B2  =  (7-7-f). 


(32) 


If  now  ud  is  assumed  to  be  given  by  a  similarity  solution  of  the  form 


=  d(0  ♦(n, ) , 
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where  tij  =  y/4,  the  equation  (32)  can  be  transformed  to 


■jP#*-'1)- 


l  a  /  J— \ 

-75^7  m 

ni  1 


Integrating  (34)  once  with  respect  to  from  n  jto  •  and  using  the  boundary 
conditions  Ud(«>)  =  7(->)  s  0,  yields 

Cni+1  ♦  +  {(J+l)  +  /  nJ1»dni]  *  ni  x  .  (35) 

5  S'q  nj  1  1 

If  now  an  eddy  viscosity  Is  assumed  to  obtain  the  turbulent  shear 
stress,  7  can  be  written 

-  vTp2  1 

t  =  2  ^  <1  ♦' »  (36) 


where  R  Is  the  Reynolds  number  based  on  the  freestream  velocity  and  length 
L.  Substituting  (36)  Into  (35)  and  rearranging  gives 


"  ~3 — T  [tii+1  *  *  {J+1  +  “  2  R1 


If  i>  Is  to  be  a  function  of  nj  only  then  It  Is  required  that 


(J+l )  +  f^q  *  (J+1 )  c5 » 


where  C5  and  Cg  are  constants.  The  right-hand  side  of  (38)  is  chosen  to  be 
(J+l)  C5  for  convenience  in  handling  both  two-dimensional  and  axisymmetric 
flows  in  a  similar  way.  If  the  eddy  viscosity  is  assumed  to  be  given  by  the 
relation 


vt/R  =  Dc"1"1  (40) 

a  solution  of  equation  (38)  and  (39)  can  be  written  down  as 

«2  ~  q!.#sSM)(W)(VD  (41) 


for  m  i  -  1  and 

62~A6lns,  q2  *  A?(ln  5)J+1(C5_1),  (42) 

for  m  =  -  1 , 

where  A4,  Ag,  Ag,  A7  are  all  constants. 


Using  the  constraints  (38)  and  (39)  in  (37)  and  differentiating  (37)  with 
respect  to  rij  once  leads  to 

Cg  4  *"  +  (Cg  +  nf )  +  (J+l)(l-C5)ni+  =  0.  (43) 


The  above  equation  can  be  transformed  to  the  well  known  confluent 
hypergeometric  differential  equation  by  changing  the  variables  to 
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Then  equation  (43)  reduces  to 

z  *J'(z)  +  (1-z)  *{{z)  -  =  0  (44b) 

where 

bx  -  C(J-l)  -  (J+l)C5]/2  .  (44c) 

A  general  solution  for  can  be  written  down  (see  Abromowitz  and  Stegun, 
1972)  in  terms  of  Kummer's  functions  as 


*  =  kjMflj.b.z)  +  k2  U(b1,l,z) 


(45) 


where  kj  and  k2  are  arbitrary  constants,  M  is  the  Kummers  function  and  U  is 
related  to  Kummers  function.  M  and  U  are  given  by 


M  (bj.l  z)  -  |  (b^,, 


(46a) 


and 


1  .  (Mr/ 

U(bl1»z)  =  rTj^T^M1’2)  1n  2  +  z  ~~ 2  C»2(b^-n)  +  l»2(l+n)]} 


(46b) 


where  i»2(bj+n)  is  known  as  the  Digamma  function,  a  logarithmic  derivative  of 
the  Gamma  function  r  (bj+n),  given  by 


♦2 (bj+n)  3  n-l+bj  +  n-2+bj  +  ••••  1+bj  +  ^  +  MM* 


(47a) 


Here 


(Mn  “  b1(b1+l)....(b1+n-l)  , 


(47b) 
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and 

♦2 (bi )  =  'Hb^  ^ln  r^bi^ •  (47c) 

For  large  z  the  behavior  of  M  and  U  are  given  by 

l  7  V1  1 

M(bl*1*z)  sTT^Te  z  El  +  0(z  X)]  (48) 

’bl  1 

and  U(b1#l,z)  ■  z  A[1  +  0(z-1)]  .  (49) 

The  above  trend  in  M  for  large  z  makes  the  velocity  profile  1(1  (see  (44a))  decay 
algebraically  at  large  for  ^  <  1  and  blow  up  for  b^  >  1.  Since  it  is 
required  that  the  solutions  decay  exponentially  to  the  freestream  values  at 
the  edge  of  the  shear  layer  kj  in  equation  (45)  must  be  set  equal  to  zero. 
Because  of  the  algebraic  behavior  of  U  for  large  z,  tends  exponentially  to 
zero  at  large  z.  The  solution  for  4>(nj)  can  be  written  down  as 

*(nx)  =  k2e'z  Ufbj.l.z),  z  *  (nf/2c|)  .  (50) 

The  behavior  of  this  solution  for  small  is  given  by 

q" (5)"  =  =  kl^_  r(b1)  *1n  nl  +  C7  +  °^nl  (51) 

Here  C7  is  another  constant.  The  above  solution  has  the  required  logarithmic 
behavior  for  small  nj  to  match  with  inner  wake  solutions  for  large  ?.  Since 
for  small  the  inner  wake  solution  gives  the  behavior  for  u 
as  u/ut  »  y  in  ni  +  CQ,  q(c)  must  be  proportional  to  ut  and  hence  must  be  a 
constant  given  by  2k *  (uTAJ#e)/k.  The  value  of  C5  therefore  turns  out  to 


22 


be  1  from  equations  (41)  and  (42).  Hence  b^  will  be  equal  to  1  for  both  two- 
dimensional  and  axisymmetric  flow  as  can  be  seen  from  (44b).  The  value  of  Cg 
can  be  arbitrarily  set  to  unity  without  loss  of  generality  because  it  only 
serves  to  determine  the  definition  of  6. 

6.  Discussion 

No  attempt  is  made  here  to  compute  the  detailed  velocity  profiles  either 
in  the  inner  wake  or  in  the  outer  wake.  Comparison  of  the  analysis  of  Alber 
with  the  two-dimensional  experimental  data  has  been  made  by  Alber  (1980)  and 
Ramaprian  et  al.  (1981).  These  comparisons  show  that,  for  the  two-dimensional 
wake  flow,  the  limiting  variations  of  the  velocity  profile  like  the  log  law  at 
the  outer  edge  of  the  Inner  wake  and  logarithmic  variation  of  the  centerline 
velocity  agree  reasonably  well  (see  figures  3  and  4)  inspite  of  the  fact  that 
Alber' s  analysis  did  not  give  a  zero  velocity  gradient  at  the  centerline.  The 
present  analysis  has  an  additional  option  of  choosing  p(c),  the  function  which 
determines  the  departure  of  the  eddy  viscosity  from  its  linear  variation 
with  c  at  the  centerline.  By  properly  choosing  p  It  is  possible  to  obtain  the 
correct  variations  of  the  velocity  profiles  near  the  centerline  in  addition  to 
the  correct  limiting  variation  for  large  c  and  large  x+.  Since  well- 
established  experimental  data  are  not  available  for  wakes  behind  an 
axisymmetric  body  and  a  yawed  plate  with  zero  pressure  gradient  fully- 
developed  turbulent  boundary  layers,  no  comparison  can  be  made  at  present, 
even  with  respect  to  the  limiting  variations  of  the  velocity  profiles  and 
their  scales. 

Since  no  predictions  have  been  made  earlier  for  the  outer  wake 
development  It  will  be  Interesting  to  examine  the  two-dimensional  wake  data  in 
the  light  of  the  analysis  made  here.  A  cursory  examination  of  the  data  of 
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Sastry  (1)  in  the  wake  of  a  flat  plate  show  that  the  total  wake  width  varies 
very  slowly  in  the  near  wake  region.  Equation  (41)  therefore  suggests  that 
(m+1)  should  be  very  nearly  zero.  If  (m+1)  is  zero  then  the  variation  of  6  is 
given  by  equation  (42)  which  shows  that  6  is  proportional  to  In  c  or  In  x  . 
Figure  5  shows  the  variation  of  the  total  width  of  the  wake  with  x+  The  wake 
width  is  defined  by  the  condition  u  =  0.995  U^.  Figure  5  shows  clearly  the 
logarithmic  variation  as  predicted  by  the  present  analysis  (equation  (42)). 
Since  (m+1)  is  very  near  zero  it  would  mean  that  the  eddy  viscosity  in  the 
outer  near  wake  region  is  constant  even  with  x  and  is  determined  by  the  outer 
layer  of  the  upstream  boundary  layer.  Departures  from  the  logarithmic 
variation  occur  when  the  maximum  defect  velocity  has  reduced  appreciably  (to 
around  20%  of  the  freestream  velocity)  and  probably  there  the  limit  considered 
in  the  present  analysis  is  no  longer  valid. 

The  agreement  of  the  predictions  with  the  experimental  data  with  respect 

to  the  variations  of  the  different  scales  and  the  velocity  profile  in  the 

overlapping  region,  suggests  that  there  Is  only  one  limiting  condition 

leads  to  the  description  of  the  near  wake  flow.  This  limiting 

u  £ 

condition,  ut/Uoo  -*•  0  as  *  *  is  the  same  as  that  used  for  an  asymptotic 
description  of  the  fully  developed  turbulent  boundary  layer.  Since  the 
maximum  defect  velocity  in  the  wake  becomes  small  at  the  end  of  the  near  wake, 
it  is  possible  to  describe  the  region  beyond  the  near  wake  by  the  linearized 
wake  equations  which  are  the  same  as  the  outer  near  wake  equations.  The 
initial  condition  for  the  solution  of  the  linearized  wake  equation  could  be 
obtained  from  the  asymptotic  near  wake  solutions  for  large  x+.  Before  the 
linearized  wake  develops  into  a  self  preserving  asymptotic  far  wake  (Narasimha 
and  Prabhu,  1972),  a  relaxation  In  the  turbulence  closure  hypothesis  from  the 
boundary- layer  type  to  the  far-wake  type  would  be  required.  The  differences 
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observed  In  the  development  of  the  scales  of  the  velocity  profile  In  this 
transition  region  (Ramaprian  et  al.  1981)  as  compared  to  the  far  wake  are 
essentially  due  to  this  relaxation  In  the  shear  stress. 

Conclusions 

Analysis  of  turbulent  near  wakes  developing  from  zero  pressure-gradient 
fully-developed  turbulent  boundary  layer  at  two-dimensional,  axisymmetric  and 
infinitely-swept  trailing  edges  have  been  considered.  It  is  shown  that  the 
same  limiting  condition  (u^/U^  ♦  0)  as  used  for  an  asymptotic  treatment  of 
the  under- determined  equations  (with  no  assumption  on  the  closure  hypothesis) 
of  a  fully-developed  turbulent  boundary  layer  for  large  Reynolds  numbers  in 
zero  pressure  gradients,  is  also  applicable  to  the  analysis  of  the  near 
wake.  It  is  argued  that  such  an  analysis  leads  to  the  governing  equations  for 
different  regions  of  the  near  wake  and  the  limiting  logarithmic  velocity 
variations  in  both  the  streamwlse  and  normal  directions  in  the  overlapping 
regions.  It  is  shown  that  similarity  solutions  can  be  obtained  for  both  the 
inner  and  outer  regions  which  tend  to  logarithmic  velocity  variations  In  the 
overlapping  region.  Comparison  with  recent  experiments  In  two-dimensional 
wakes  behind  a  flat  plate  shows  that  the  most  relevant  solutions  for  the  outer 
wake  corresponds  to  that  In  which  the  square  of  the  total  wake  width  varies 
slowly  like  log  x  In  the  near  wake.  Finally,  It  Is  argued  that  the  flow  In 
the  region  beyond  the  near  wake  can  be  obtained  by  using  linearized  far-wake 
equations  with  the  initial  conditions  for  the  velocity  and  shear  stress 
obtained  from  the  near-wake  solution. 
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